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Electromagnetic properties provide information about the structure of strongly interacting sys- 
tems and allow for independent tests of hadronic models. The radiative decay of the A* (1700) is 
studied, which appears dynamically generated in a coupled channel approach from the rescattering 
of the (3/2 + ) decuplet of baryons with the (0~) octet of pseudoscalar mesons. The radiative decay is 
predicted from the well-known couplings of the photon to the mesons and hadrons which constitute 
this resonance in the dynamical picture. 
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I. INTRODUCTION 



PACS numbers: 24.10.Eq, 25.20.Lj, 11.30.Rd 

o 
o 

The unitary extensions of chiral perturbation theory (UxPT) have brought new light to the meson-baryon in- 
teraction, showing that some well-known resonances qualify as being dynamically generated. In this picture the 
Bethe-Salpeter resummation of elementary interactions, derived from chiral Lagrangians, guarantees unitarity and 
leads at the same time to genuine non-perturbative phenomena such as poles of the scattering amplitude in the com- 
plex plane of the invariant scattering energy y/s, which can be identified with resonances. Coupled channel dynamics 
' plays an essential role in this scheme, with the chiral Lagrangians providing the corresponding transitions of the 
, multiplets; even physically closed channels contribute as intermediate virtual states. 

After earlier studies in this direction explaining the A(1405) and the A''* (1535) as meson-baryon (MB) quasibound 
states [1-6] from the interaction of the meson octet of the pion (M) with the baryon octet of the nucleon (B), new 
efforts have been undertaken [7, 8] to investigate the low lying 3/2 _ baryonic resonances which decay in s-wave into 
CP mesons (M) and 3/2 + baryons (B*) of the decuplet. The leading interaction of these hadrons, given by the 
isovector term from Ref. [9], is unitarized by the use of the Bethe-Salpeter equation (BSE) in the on-shcll reduction 
scheme from [8] which allows for a factorization of vertices and the intermediate loop function, thus reducing the BSE 
to an algebraic matrix equation in coupled channels. The unitarized amplitude develops poles in different isospin and 
strangeness channels in the complex plane of y/s, which have been identified with resonances from the PDB [10] such 
^ i as the A*(1520), E*(1670), A*(1700), etc. [8]. 

In isospin 3/2, strangeness 0, which are the quantum numbers of the A* (1700), the coupled channels in s-wave are 
given by A(1232)tt, E*(1385)AT, and A(1232)t/. The A*(1700), together with a series of other production mechanisms, 
has been included in Ref. [11] in the study of the jp — > ir°rip and jp — ► 7r°A'°I] + photoproduction reactions, currently 
measured at ELSA/Bonn. In the detailed work of [11] the A*(1700), together with its strong couplings to A(1232)?7 
and E*(1385)A', turned out to provide the dominant contribution. The branching ratios into these two channels are 
predicted from the scheme of dynamical generation and differ from a simple ST/ (3) extrapolation of the A-7T channel 
by up to a factor of 30 [12]. 

The predictions for both reactions are in good agreement with preliminary data [13]. Recently, new measurements 
at low photon energies have been published [14] which also agree well with [11]. This has motivated another study 
[12] of altogether nine additional pion- and photon-induced reactions. From considerations of quantum numbers and 
the experimentally established s-wave dominance of the E* production close to threshold, the A* (1700) channel is 
expected to play a major, in some reactions dominant, role. Indeed, good global agreement has been found for the 
studied reactions that span nearly two orders of magnitude in their respective cross sections. 

Thus, evidence from quite different experiments has been accumulated that the strong A* (1700) — > E*A, Ar) cou- 
plings, predicted by the coupled channel model, are realistic. This gives support to the scheme of dynamical generation 
of this resonance. However, in all the photon-induced reactions from [11, 12] the initial -fp — > A* (1700) transition 
has been taken from the experimental [10] hclicity amplitudes Ai/ 2 and ^.3/2 [15]: In Ref. [16] the electromagnetic 
form factors G[, G' 2 , and G' 3 , which appear in the scalar and vector part of the jp — > A* (1700) transition, have been 
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expressed in terms of the experimentally known Ax/2 an d A 3 / 2 [10]; this provides the transition on which we rely in 
all the photoproduction reactions via A* (1700) in [11, 12]. 

Such a semi-phenomenological ansatz is well justified: the photon coupling and the width of the A* (1700) is taken 
from phenomenology, whereas the strong decays of the A* (1700) into hadronic channels are predictions from the 
unitary coupled channel model; the strengths of these strong transitions are responsible for the good agreement with 
experiment found in [11, 12]. It is, however, straightforward to improve at this point, and this is the aim of this study. 

Electromagnetic properties provide additional information about the structure of strongly interacting systems and 
allow for an independent test of hadronic models, in this case the hypothesis that the A* (1700) is dynamically 
generated. A virtue of the present model is that one can make predictions for the radiative decay, or equivalently, the 
inverse process of photoproduction; the components of the A* (1700) in the meson-baryon base are all what matters, 
together with the well-known coupling of the photon to these constituents. A similar study has been carried out for 
the radiative decay of the A* (1520) in Ref. [17] that has been described recently as a dynamically generated resonance 
[18]. For the A* (1520) — > S°7 decay, where the dominant channels add up, indeed good agreement with experiment 
has been found. The present study is carried out along the lines of [17] but several modifications will be necessary: 
the (ttN) channel in d-wave plays an important role and is implemented in the coupled channel scheme. Second, a 
fully gauge invariant phototransition amplitude for the s-wave channels is derived that includes also couplings of the 
photon to the B* baryons. 



II. THE MODEL FOR THE RADIATIVE A* (1700) DECAY 



In Sec. II A the coupled channel model from Ref. [8] is revised and extended to the inclusion of the ttN channel 
in rf-wavc, (irN)d- In Sec. II B the model for the phototransition amplitude A* (1700) — > 7 A is derived: The photon 
interacts with the dynamically generated resonance via a one-loop intermediate state that is given by all coupled 
channels which constitute the resonance. 



A. The (irN)d channel in the unitary coupled channel approach 

In the unitarized model from Ref. [8] the A* (1700) resonance appears as a quasi-bound state in the coupled channels 
Air, T,*K, and Arj. The attraction in these channels leads to the formation of a pole in the D 33 channel which has 
been identified with the A* (1700) resonance. In Ref. [8] the chiral interaction from [9] is adapted in a nonrelativistic 
reduction providing isovector MB* — > MB* transitions where M (B*) stands for the octet of 0~ pscudoscalar mesons 
(3/2 + decuplet baryons). This interaction is unitarized by the use of the Bethe-Salpeter equation (BSE) 



T = (1 - VG^V 



(1) 



which turns out to be a matrix equation in coupled channels, factorized to an algebraic equation according to the 
on-shell reduction scheme of [8, 19]. In the recent work of [17] this scheme is compared in detail to other possible 
treatments of the BSE. The function G is a diagonal matrix with the MB* loop functions Gmb* of the channels i 
which are regularized in dimensional regularization with one subtraction constant a. As it will appear in a different 
context in the phototransition, the function Gmb* has been re-derived with the result 
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where m (M) is the meson (decuplet baryon) mass, q 2 = s is the invariant scattering energy, and 



"(A*) =7 log (4tt) 



(3) 



with the Eucledian integration over £e and e = 4 — d. The cm. energy y/s in Eq. (2) and all the following expressions 
of this study has to be taken at the physical sheet, i.e., yfs — > \fs + ie. The only exception is the pole search in the 
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second Ricmann sheet discussed at the end of this subsection. Values of the rcgularization scale of /i = 700 McV and 
a = —2 are natural, as argued in [8]. The cm. momentum function Q and f\ in Eq. (2) are given by 



V( 8 -(M + m)2)(s-(M-m)2) 

QWS) = ^ , 

= log (a - (M 2 - m 2 ) + 2yf~s Q{yfs)) + log (s + (A/ 2 - m 2 ) + 2yfs Q{y/s)) 

- log (-s + (M 2 - m 2 ) + 2y/7s Q(V*)) - log (-s - (M 2 - m 2 ) + 2^~s Q{^fs)) . (4) 

The loop function from Eq. (2) has a real part, which implies a major difference of the present approach compared 
to unitarizations with the K -matrix. The real parts of the d, together with the attractive kernel V in the isospin 
3/2, strangeness channel, provide enough strength for the formation of a pole in the complex plane of the invariant 
scattering energy s/s which is identified with the A*(1700). 

However, additional channels will also couple to the dynamically generated resonance, changing in general its 
position and branching ratios, as these new channels can rescatter as well. In this study, the (nN)^ channel is 
included in the analysis, because this is the lightest channel that can couple to the A* (1700) and precise information 
of the ttN — > ttN transition in the D33 channel exists from the partial wave analysis (PWA) of Ref. [20]. The {pN) s 
channel has been found important [10, 21], but for the radiative decay its influence is expected to be moderate as 
discussed below. 

In order to include the (irN)d channel in the coupled channel model, one has to determine the (TrN)d — ► (MB*) S 
transitions, where MB* stands for the channels An, T,*K, and A77 from [8]. There is no experimental information 
on these transitions. From the theoretical side, there is no information either due to the large number of low energy 
constants in the <i-wave to s-wave transition. Thus, the coupling strengths have to be introduced as free parameters, 
called (3. For the inclusion of d-wave potentials, we follow the lines of Ref. [18, 22] where it has been shown that 
the d-wave transitions can be factorized on-shell in the same way as the s-wave transitions; as a consequence, the 
mcson-baryon d-wave loop function is the same as the s-wave loop function from Eq. (2) . We also allow for a d-wave 
to d-wave transition of the (nN)d channel. 

With the channel ordering i = 1 • ■ • 4 for An, T,*K, Ary, (nN)d the interaction kernel is given by 



V = 



2 



Qi<N r P(KN) d ->A v 
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(5) 



where we have multiplied some elements with r = 1/(4 • 93 2 ■ 1700) MeV~ 3 = 1.7 • 10 -8 MeV~ 3 in order to obtain 
dimcnsionless transition strengths (3 of the order of one. In Eq. (5) one can recognize the chiral isovector transitions 
with (k° + k'°) from Ref. [8] where k° = (s - M 2 + m 2 )/(2^s) is the meson energy and f„ = 93 McV. In Eq. (5), 
QirN is the on-shell cm. momentum of the ttN system and the (3 are the s-wave to d-wave transition strengths. See 
also Refs. [18, 22] where the analytic form of the transitions Va (i = 1, • • • ,4) has been derived. 

Although a natural value for the subtraction constants is given by a = —2 [8], it is a common procedure [5, 18] to 
absorb higher order effects in small variations around this value. Thus, as these higher order effects are undoubtly 
present, we allow for variations of the a of the four channels. Together with the transition strengths f3, the set of free 
parameters is fitted to the singlc-energy-bin solution of the PWA of Ref. [20] . Note that there is a conversion factor 
according to 



^--Vs^gjW!) (6) 

in order to express the solution T of the Bethe-Salpeter equation (1) in terms of the dimensionless amplitude T(y/s) 
from [20]. In Eq. (6), Mi {Qi) is the baryon mass (cm. momentum) of channel i; in the present case, i = j = 4. 

For energies above 1.7 GeV a theoretical error of 0.08 has been added to the error bars from [20] because additional 
channels such as pN start to open and one can not expect good agreement much beyond the position of the A* (1700). 
The resulting amplitudes of four different fits are plotted in Fig. 1 with the parameter values displayed in Tab. I. 

In fit 1, all parameters have been left free. The values for the (3 are small: For y/s ~ Ma*, they lead to values in 
the kernel (5) around one order of magnitude smaller than the chiral interactions between channels 1 to 3; the size of 
(3(-nN) d ->(nN) d corresponds to a value sometimes even two orders of magnitude smaller. These parameter values reflect 
the fact that the A* (1700) couples only weakly to itN and that the itN interaction in D33 is weak in general. 
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FIG. 1: Fit results: Fit 1: All parameters free. Fit 2: without ttN — > nN d-wave transition kernel. Fit 3: As fit 2, but all 
subtraction constants for the s-wave loops chosen to be equal. Fit 3': As fit 3 but different minimum in \ ■ The error bars 
show the single-energy-bin solution T44 from Ref. [20]. 



TABLE I: Parameter values of the fits to (nN)d'- the on are subtraction constants, the fli are (irN)d —> (B* M) 3 and (nN)d — > 
(irN)d transition strengths. Note the sign changes of the Pi between fit 3 and fit 3'. 
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Fit 1 


— 1.96 


-1.28 


-0.87 


-1.00 


2.19 


-1.15 


-0.07 


63.5 


Fit 2 


— 1.25 


-1.29 


-0.66 


-1.96 


3.30 


-0.37 


0.31 





Fit 3 


-1.21 


-1.21 


-1.21 


-1.00 


3.25 


-0.85 


-0.54 





Fit 3' 


-1.22 


-1.22 


-1.22 


-1.00 


-3.19 


0.92 


0.42 






Thus, the (irN)d — > (TrN)d transition strength can be set to zero which is done for the fits 2, 3, and 3'. One can 
even choose all subtraction constants of the s-wave channels to be equal, which is done in fit 3, and still obtain a 
sufficiently good result as Fig. 1 shows. However, for fit three, there is another minimum in x 2 , almost as good as 
the best one found. This fit is called fit 3'. As Tab. I and Fig. 1 show, one obtains an almost identical amplitude 
with a set of /3's with opposite sign as compared to fit 3 (see explanation below). 

For the different solutions 1 to 3', the coupling strengths of the resonance to the different channels can be obtained 
by expanding the amplitude around the resonance position in a Laurent series (see also Sec. HE). The residues give 
the coupling strengths which are uniquely determined up to a global sign which we fix by demanding the real part of 
the coupling to An to be positive. In Tab. II the resulting couplings are displayed. The values in brackets quote the 
values of the original model from Ref. [8] without the inclusion of (irN)d- Compared to these values, all couplings 



TABLE II: Position and couplings of the A* (1700). The values in brackets show the original results from [8] without 

the inclusion of the (7rA) d -channel. The PDB [10] quotes a value of = (1620 - 1680) - i (160 - 240) MeV and couplings 
corresponding to IpattI = 1-57 ± 0.3, \g( W N) d I = 0.94 ± 0.2. Note the sign change of g(TrN) d between fit 3 and fit 3'. 
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Fit 1 


1707 - i 160 


1.09 - i 0.92 


1.4 


3.57 + i 1.91 


4.0 


-1.98 


- i 1.68 


2.6 


-0.84 - 


i0.05 


Fit 2 


1692 - i 166 


0.62 -i 1.03 


1.2 


3.44 + i 2.28 


4.1 


-1.89 


- i 1.78 


2.6 


-0.77 4 


iO.OO 


Fit 3 


1697- i 214 


0.68 -i 1.07 


1.3 


3.01 +i 1.95 


3.6 


-2.27 


— a 1.89 


3.0 


-0.89 4 


i0.15 


Fit 3' 


1698 — i 216 


0.68 -i 1.07 


1.3 


3.02 +i 1.95 


3.6 
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- i 1.89 


3.0 


40.92 - 
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increase slightly in strength and have different phases. The largest change is observed for I^attI which has increased 
by 40 % for fit 1; the new value is well inside the range quoted by the PDB [10] that corresponds to |<7att| = 1.57±0.3. 
Also, the d-wave coupling to irN coincides well with \g(nN) d I = 0.94 ± 0.2 [10] . The main properties of the resonance 
are conserved; in particular, the absolute values \g\ do not change much. This is relevant with respect to previous 
studies [11, 12] where the model from [8] has been used for the couplings of the A*(1700) to Ar) and T,*K. 

Some explicit corrections to the results from [11, 12] are given, based on the values from fit 1 which is the preferred 
one as discussed below. The cross sections of the pion induced processes from [12] with the T,*K final state will 
change by a factor (| — 0.84 — i 0.05|/0.94) 2 (4.0/3.4) 2 = 1.1, i.e., they stay practically the same. In the pion induced 
reaction with the A77 final state the factor is again 1.1. Anticipating the result for fit 1 from Tab. IV, that the 
radiative coupling is close to the experimental one (which is the one used in [11, 12]), the cross section of the reaction 
7P — > A'°7r + A increases by a factor of (4.0/3.4) 2 = 1.4 which improves the agreement with data, see Fig. 5 of [12] 
and also the recent measurements in Ref. [23]. For — > ir°K°'E + , which is also studied in [11, 12], and where the 
photoproduction via the A*(1700) dominates [11], the factor is also close to 1.4. For the cross section of the 7p — > TT°T]p 
reaction, the transition via the A* (1700) is only one of the reactions, and one has to re-evaluate the coherent sum of 
all processes of the model of [11], with the new values from Tab. II. The numerical results for this reaction [11] have 
been updated in Fig. 5 of [12]. With the values of fit 1 from Tab. II, the cross section stays practically the same as 
in Fig. 5 of [12] (10 % decrease). 

Interestingly, the sign of the coupling to (ttN)^ is reversed in fit 3' compared to fit 3. This behavior has been 
noted before for the parameter values in Tab. I. The reason for the difference between fit 3 and 3' is that the 
coupling of the irN channel to the A* (1700) is small: Although (irN)d is contained to all orders in the rescattering 
scheme, higher orders are smaller. Then, the (by far) dominant order in the fit to (irN)d — > (irN)d is gfnN) d an< ^ the 
relative sign to the other couplings 5Att, 5s*k, and <?a?j is difficult to fix. As the influence of the (irN)d channel in 
the rescattering scheme is moderate, the main differences to the original model arc caused by other parameters: the 

important parameters, which have fine-tuned the original model from [8], are the subtraction constants oa^, cce*k, 

1/2 

and a Art- Their variation brings the pole down from s olc = 1827 — i 108 McV [8] to the values quoted in Tab. II. 

In Sec. Ill results for the radiative decay for all four fits are given, in order to obtain an idea of the systematic 
theoretical uncertainties. The fit 1 is preferred, though, because in the reduction to less free parameters, as it is the 
case for the fits 2, 3, and 3', the remaining free parameters have to absorb effects such as direct (irN)d — * {nN)d 
transitions; results might become distorted, despite the fact that the fits appear to be good in Fig. 1. The larger 
space of free parameters also helps to fix the ambiguities found in fit 3 and 3': For fit 1 and 2, no alternative solutions 
with a reversed sign for g( n N) d have been found, and, thus, the sign is fixed. 

I have also performed a search for poles in the second Riemann sheet. For the different fits, the pole positions 
are given in Tab. II. The positions are in agreement with the values given by the PDB [10]. Note that one of the 
virtues of the coupled channel analysis is that a separation of background and resonance part of the amplitude is not 
necessary; thus, the position of the pole does not suffer from ambiguities from this separation process required by 
other analyses. 

Nevertheless, some theoretical uncertainties are present in the model from the omission of other channels such as 
pN in s-wave or even pN, An and KY, in d-wave, all of which have been reported in the PDB [10]. Although the pN 
channel is closed at the position of the A* (1700), it contributes through the real part of the pN loop function in the 
rescattering scheme, and through the finite width of the p even to the imaginary part. 

However, in the calculation of the radiative decay, which is the aim of this study, no large contributions are expected 
from these heavy channels. In the study of the radiative decay of the A* (1520) [17] we have seen that contributions 
to the radiative decay width from the heavy channels are systematically suppressed, as also discussed in Sec. III. 

B. The phototransition amplitude 

The only known radiative decay of the A*(1700), which is also the one of relevance in Rcfs. [11, 12], is into 7A 
and we concentrate on this channel. The coupled channel model for the A* (1700) has the virtue that the radiative 
decay can be calculated in a parameter-free and well-known way through the coupling of the photon to the particles 
which constitute the resonance. The dominant photon couplings to the coupled channels Att, Y?K , A77, and (irN)d 
are displayed in Fig. 2 . We have chosen here the ^3 = +1/2 or charge C = +1 state of the A*(1700). The set of 
diagrams (1) to (9) is gauge invariant and finite as shown in Sec. II C. 

There is also a loop diagram where the photon couples to the A + in an intermediate A + r/ state. This contribution 
is doubly suppressed: First, the diagrams with ^B*B* couplings arc much smaller than the other ones (see Sec. II C), 
and, second, the At? state is heavy which renders this contribution even smaller. 

Apart from the photon coupling to the An and T,*K channels, Fig. 2 shows also the coupling to the irN d-w&ve 
channel in diagram (10). Note that there is no Kroll-Ruderman term (which has irN in s-wave in the 7-BMB-vertex) 
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FIG. 2: Mechanisms for the A* (1700) decay in s- and d-wave loops. The shaded circles represent the A* (1700) in the transition 
B* M — + An or -K + n — > Ai. The diagrams in the left column are referred to as "meson pole loops". Diagrams (2), (5), and 
(8) have a 7BMB* transition and are referred to as " Kroll-Ruderman loops" (diagrams (3), (6), (9): "baryon pole loops"). 
Diagram (10) shows the process with a d-wave coupling of the (irN)d intermediate state to the A* (1700). 



as the vertex on the right hand side of the irN loop is in <i-wave, and the term vanishes in the integration over the 
loop momentum. The photon coupling to the <i-wave loop is evaluated in Sec. II D. 

As for additional photon couplings, as e.g. to the external proton, to vertices of the rescattering, or to intermediate 
meson-baryon loops, apart from the ones considered, these processes are present [24] in general but negligible as 
discussed at the end of Sec. II C. 

The AIBB* vertices appearing in Fig. 2 are provided by the Lagrangian from Ref. [25], with the relevant parts for 
the present reaction given by 

C = C (T^B + BA fl T f *) 

(!,■■■ ,3 l,-,3 \ 

E t abc T ade u,A b ^ B<+ £ e abc B e c <»T ad e U »\ 
a,b,c,d,e a,b,c,d,e J 

(7) 

where B is the standard 3x3 matrix of the l/2 + baryon fields, T is the field of the 3/2 + decuplet baryons with the 
definitions and phase conventions of [25], and A M is proportional to the axial current [25]. In the present situation, A 11 
is expanded up to one meson field. From minimal coupling of the photon, the ^BMB* Kroll-Rudcrman contact term 
is straightforward deduced from Eq. (7). The constant C is fitted to phcnomcnologically known branching ratios. 
Explicit expressions for the required Feynman rules can be found, e.g., in [11, 12, 17]. 



7 



Let us start with the evaluation of diagram (1) in Fig. 2. The amplitude is given by 

m V o M 

x (-*) I 77TTX1 W^TJ-o --o K -7 T~vT~ 5 ~ (p - kY{2p - k) v (8) 

v ' J (2tt) 4 (q - p) 2 - M 2 + it p 2 - m 2 + it (p - k) 2 - m 2 + le y ' y ' w 

where m(M) is the pion (A) mass, f A ^N = m n /(V2U)C a ^n-k = 2.13 is the AirN coupling strength, e 2 = Ait/131 
is the electric charge, is the spin 1/2 — > 3/2 transition operator which we approximate by = (0,S>), and e v 
is the polarization of the photon which in Coulomb gauge is given by t v = (0,e). The shaded circle in diagram 
(1) represents the T-matrix element T^-^a-k of the unitary coupled channel scheme in which the A* (1700) appears 
dynamically generated. In Sec. HE it will be matched to the results from Sec. II A. Note the simplified structure of 
the A-propagator in Eq. (8). This is the same simplification as made in the model of dynamical generation [8] of the 
A*(1700), see Eq. (2). 

In order to ensure gauge invariance, the loop in Eq. (8) can be evaluated using a calculation technique from Rcfs. 
[26, 27]. In Sec. II C we will compare this scheme to a straightforward calculation of the loops of Fig. 2. The general 
structure of the loop function, or phototransition amplitude, is given by 

T H» = ag ^ + bq»q u + cq»k v + dk^q" + ek»k v . (9) 

with the momenta q and k as defined in diagram (1). The terms with c and e do not contribute once contracted with 
e" according to Eq. (8) and using the transversality of the photon, efc = 0. The terms with b and d do not contribute 
as eq = in the cm. frame where |q| = and using the fact that e° = 0. Thus, the only term that will not vanish in 
Eq. (9) is (ag^). 

It can be shown that the sets of diagrams (1) to (6) and (7) to (9) of Fig. 2 are gauge invariant. Contracting T M1/ 
with the photon momentum k v and using the Ward identity k y = 0, leads to the condition a + dkq = 0. Note 
that diagram (2) of Fig. 2 contributes only to the term with a, whereas diagram (1) contributes both to a and d. 
However, evaluating d (and from this, a through the condition a + dkq = 0) has the advantage that the loop integral 
is finite whereas both diagram (1) and (2) are logarithmically divergent. Using Feynman parameters and keeping only 
the terms proportional to k^q v , the second line of Eq. (8) becomes 

„_ 4Mk»q" f V x(z-l) 

q (4tt 2 ) J J x[(x - l)q 2 + z(q 2 - M 2 ) + M 2 ] + (1 - x)m 2 ' 1 J 

o o 

where we have written the product 2qk = 2q°k° = q 2 — M 2 in the cm. system where |q| = and M e is the mass of 
the external baryon, in this case a proton. Note that fc° = |k| = l/(2 v / i)(s — M 2 ) where yfs = q° which we will use 
several times in the following. From Eq. (10) we calculate the term a through the condition a + dkq = 0, 



G a ; = a = —dkq 
2M 



kJ~s \(M 2 - m 2 ) 2 - 2M 2 M 2 ] - m 2 M* M 2 

~ lQ g — 



2A/ 2 2 Mi s TTL 



(M 2 - M 2 + m 2 ) (2M 2 - s) M 2 - M 2 -m 2 + 2k y T s 

iMjfcTi Q(Mc) h{Mc) + Ik's Q{VS) h{Vs) 

- 2 r T . f-M 2 + M 2 + m 2 -2M C Q{M C )\ ( -M 2 + M 2 + m 2 + 2M B Q(M e ) 

-Ll2 Z 5 + 1^2 



2k y/s 



2m 2 J \ 2m 2 



:-M 2 + s + m 2 -2^fsQ{^fs)\ ( -M 2 + s + m 2 + 2^/s Q(s/s) \ 

Ll2 1 2^ r 2 \ ^ ) 



(11) 



which is gauge invariant by construction. In Eq. (11), Li2 is the dilogarithm and Q and /i are given in Eq. (4). 
Having calculated a from d, the loop function Gg ; corresponds to the meson pole diagram (1) from Fig. 2 plus the 
Kroll-Ruderman term from diagram (2). 

Furthermore, the photon can also couple directly to the baryon as displayed in diagram (3). Only the non-magnetic 
part of the <yB*B* coupling (convection term), given in [16], is considered. This term has the same structure and sign 
as the 7MM vertex. For the magnetic part, see the discussion at the end of Sec. II C. Diagram (3) also contributes 
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to the term dk^q v in Eq. (9). The contribution to d, let it be d 11 , leads to an extra modification of the term a, 
G" = a 11 = -d n kq 



2M 
(4^r)2 

/™2 



Ml + M 2 - m 2 {M 2 - m 2 f 



2M? 

M 2 ) (2M 2 -s)-sM 2 



2MU 



, M 2 
log — 



4M| ft y/s 



Q(M C ) A (M e ) 



4fcs 



M 1 



2ky/s 
-Li 2 



Li 2 



-m 2 + M 2 + M 2 - 2M C Q(M e ) 



2M 2 



—m + 



.s + M 2 -2 v / iQ(y / i) \ 



2M 2 



Li 2 



Li9 



— m + s 



Q(V5)/i(Vi) 

M e 2 + M 2 + 2M e Q(M e ) 



2Af 2 

M 2 + 2VsQ(Vs) \ 



2M 2 



(12) 



In order to determine the effective coupling of the photon to the A* (1700) we construct isospin amplitudes from the 
diagrams given in Fig. 2. The isospin states of Air and T,*K in (/ = 3/2, 1% = 1/2) are given by 



|A7r,I = 3/2,J 3 = l/2) = W-|A 



,++, 



\X*K,I = 3/2,/ 3 = l/2) = \/-\X*+K ) 



— IA+7T ) 

15 1 7 
Z | E *0#+) 



'— |A tt + \ 
15 



(13) 



with the phase convention \ir + ) = — 11, 1). 

With the loop functions from Eqs. (11) and (12) and standard Feynman rules [11, 12, 17] we can calculate the 
isospin amplitudes for the sum of all A7r-loops and the £*AT-loops according to Eq. (13) with the result 



-it ■ e) 



(7=3/2,13=1/2) 

7p— >A7T— *A7T 



Ht . €) ('=3/V3=l/2) 



ViO 


eSt 




3 






1 


eSt 




3V2~ 







( G g.i- 



NK 



'r 1 4- r 11 ) 



^7P->£»Jf->-A7r ;a.y \ m = mK t M=M S , , M C =M N 

where we have indicated which masses m, M, M e have to be used in the loop functions. In Eq. (14) 



C^NK = 1-508 : 



6(D + F) 



(14) 



(15) 



The strength C^^mk f° r t ne ^* decay into the physically closed channel NK has been determined from from a SU (6) 
quark model [28] in the same way as in Rcfs. [11, 12]: the SU(6) spin-flavor symmetry connects the ttNN coupling 
strength to the 7riVA strength, and then SU(3) symmetry is used to connect the 7riVA transition with KNU*. The 
use of SU(6) symmetry allows to express C^^nR m terms of D and F. 

In the decuplet, the SU(3) symmetry is broken. This can be taken into account phenomenologically by allowing 
for different C in the Lagrangian (7). For the open channels of the S* decay modes one obtains Ce*->Ett = 1-64 and 
Cs»^A7r = 1-71 from fitting to the partial decay widths into these channels. The constant C-^-^nk = 1-508 from Eq. 
(15) is close to these values (compare to Ca->Ntt = IattN — 2-13). 



C. Gauge Invar iance 



The construction of the gauge invariant amplitude in the last section can be compared to a straightforward calcu- 
lation of the diagrams (1) to (9) in Fig. 2. In this section we show that both ways give identical results; at the end 
of this section we discuss further issues related to gauge invariance. 

In Fig. 2 there are three types of loops: the Kroll-Ruderman structure, the meson pole term, and the baryon 
pole term. All loops are logarithmically divergent and we calculate in dimensional rcgularization for the sake of 
conservation of gauge invariance. The Kroll-Rudcrmann loop function is identical to the common meson-baryon loop 
function from Eq. (2), 



(16) 
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TABLE III: Coefficients At for the diagrams (1) to (9) from Fig. 2 with the amplitude given in Eq. (20). The lower row 
shows how the infinities 2/e (e = 4 — d, see Eq. (3) for e — > 0) scale for each diagram. Once multiplied with the corresponding 
Clebsch-Gordan coefficients (CG) according to Eq. (13), the sum over the infinities cancels, E, (CG) t rt (-) = 0. 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 

I / 1 /"i /l/O /2/o /-~» o/o / 1 /"» / 1 /"f / 2 

-K!) 1(1) 7f(!) -#(!)V3(|) -V3(|)-I(|) 1(f) 



With the momenta assigned as in diagram (1) of Fig. 2, we define the meson pole loop function by 
■\ f di P 2M 1 1 i.y*fo 

— » G 1 MM 

g^M (2/e - 7 + log(47r) + log ^ 2 ) 2M 



. dx dz log (x[(x - 1)<? 2 + 2zqk + M 2 } + (1 - a;)m 2 ) 
{Any {Any J J 



^(^2^K M ) ( 17 ) 



with 



1/ , to 2 \ (A/ 2 - m 2 ) 2 + sM 2 , Af 2 
i?(m,M) = - (l-a-log— ) T -i 73 c - log — 



2 V M 2 / 4sM| b m : 

2 r /-Af 2 + Af o 2 + TO 2 -2Af Q(M c )\ . / -A/ 2 + A/ 2 + m 2 + 2M e Q(M e ) 



2fcVi 
—Lis 



Li 2 : — e —r— n " v w + Li 



2m 2 / \ 2m 2 

-M 2 + s + m 2 -yiQ(vj) \ / -M 2 + g + m 2 + 2yfiQ(y/s) \ 

2m 2 J 12 I 2m 2 J 



(18) 



and Q, f\ from Eq. (4). The arrow in Eq. (17) indicates that we only keep the terms proportional to g^ v because 
all other possible structures from Eq. (9) do not contribute as commented following Eq. (9). Similarly, and with 
the assignment of momenta as in diagram (3) of Fig. 2, the loop function where the photon couples directly to the 
baryon, is given by 

J" GyB*B* = gT 7^2 R(M, m). (19) 

Note that the convection part of the 7-B*i?* coupling (the non-magnetic part) is of the same structure and sign as 
the jMM coupling [16]. 

The next step is to express the amplitudes of the diagrams in Fig. 2 in terms of these three loop functions. Using 
standard Feynman rules [11, 12, 17], we obtain for the diagrams (1) to (9) 

(-it ■ e)(j) = AiegB'MB & • cTb'M^a-k (20) 

where Qb*mb = Iattn/ 711 ^ or 9b*mb = Cs*^7Vi?/( 2 .Ar) f rom E q. (15) depending on whether (B*M) = (Air) or 
(Y,*K). The coefficients A{ for the diagrams (i) for i = 1 to 9 are given in Tab. III. 



The sum over all diagrams results in 



(21) 



1 „•+ ^\( /=3 / 2 ^3=i/2) oS! t , ^s*->Arg m j.r> <r> \ t 
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The infinities of the nine diagrams, i.e., the terms with 2/e from the Eqs. (16, 17, 19) scale as shown in the lower 
row of Tab. III. In order to construct the isospin states according to Eq. (13) each infinity is multiplied with the 
corresponding Clebsch-Gordan coefficient; as a result, the sums over all infinities cancel for the Aw loops and also 
for the Y,*K loops. Note that the phase convention \ir + ) = — 11,1) is crucial at this point. The Ward identity is 
working. In other words, gauge invariance renders the phototransition amplitude finite and leads to a parameter-free 
expression. Comparing Eqs. (21) and (14), obviously 

GjMM + GjBMB* + GyB'B* = G g L + G g 7 ; (22) 

which can also be seen by comparing the explicit expressions given in Eqs. (11, 12, 16, 17, 19). In other words, 
the scheme from Eq. (9), which allows for the construction of a gauge invariant amplitude through the condition 
a + dkq = 0, leads to the same result as a straightforward calculation of the amplitude, in which the infinities cancel 
systematically. However, this is only the case if all contributions to d are taken into account, in the present case from 
the meson pole term plus the baryon pole term (G g ; and G^ ). 

In the rest of this section, further issues of gauge invariance are discussed such as a comparison to cut-off schemes, 
the role of magnetic couplings, and gauge invariance in the context of the rescattering scheme. 

In several recent studies [11, 17] the occurring photon loops have been regularized with a cut-off. As we have now 
a gauge invariant, parameter- free scheme at hand, we would like to compare both methods numerically. As a first 
test, the scheme has been implemented in the calculation of the radiative decay width of the A* (1520) from Ref. 
[17]. This means a gauge invariant evaluation of the s-wave loops from Fig. 2 of [17] formed by 7rE* and KE.*, 
plus additional diagrams with 7£*£* and 7S*S* couplings in analogy to the diagrams in Fig. 2. In practice, the 

re-calculation only requires the replacement of the terms (^G L + | Gj\ from Eq. (41) of [17] by (G g ; + ) or 

Gjmm + GyBMB' + GyB'B' from Eq. (22). The final result from Ref. [17] for the radiative decay A*(1520) — > 7E 
changes from T = 60 keV [17] to T — 61 kcV. Thus, the approximations made in [17] and the violation of gauge 
invariance are well under control. Note that the additional diagram with a 7s* £* coupling cancels for the second 
radiative decay studied in [17], A* (1520) — > 7A, in the same way as the 7rE* meson pole and Kroll-Rudcrman terms. 
This is a consequence of gauge invariance but can be also seen directly by noting that the 7MM interaction and the 
convection term of the 7 interaction have the same structure and sign [16]. Thus, the A* (1520) — > 7A radiative 

width stays as small as already found in [17]. 

In the cut-off scheme, the meson pole loop is defined as 



G (cut) = i 



A q 2 -(q-k) 2 /|k| 2 1 Mi 1 



(2tt) 4 (q - k) 2 - m 2 + ie q 2 - m 2 + ie £;(q) P° - q° - E t (q) + ie 
A 1 

dqq 2 f, q 2 (l~x 2 ) 1 1 Mi 



dx 



(2w) J luiiU)': k + uii + u)'. k — uji — u[ + ie Ei(q) 

-1 



(ui+ujtf + (u i + u>' i )(E i (q)-^)+ku>' i . (23) 



i/s-Wi - E t (q) + ie ^/s - k - - E^q) + ie 

For consistency, the notation is as in [17]; x is the cosinus of the angle between q and k with k the momentum of 
the real photon (|k| = k); rrii is the meson mass, P° = y/s, and u>i, are the energies of the mesons at momentum q 
and q — k, respectively; Ei the energy of the baryon. Note that Eq. (23) is slightly different from the corresponding 
expression in [17], called G there, as one kinematical approximation made in [17], the angle average over the term 
1 — x 2 , is not sufficient in the present case, because the momenta of the decay products are higher. One obtains G 
from [17] by substituting (1 - x 2 ) -> 1 in Eq. (23). Note that G (cut) corresponds to 2/3 of G. 

The meson pole term G( cut ) is accompanied by the corresponding Kroll-Ruderman term with cut-off, called G^ cnt ^ , 
where explicit expressions can be found, e.g., in Eq. (22) of [17]. The cut-off can be determined by requiring the real 
part of the Kroll-Ruderman loop function to be equal in both dimensional regularization and cut-off scheme, at the 
energy of the resonance. For An in the loop, this leads to a cut-off of A = 881 MeV. 

In Fig. 3 the cut-off loops for An (G' cut ' + G^ cut - ) ), are shown as the dotted line. The gauge invariant function 
Gg ; from Eq. (11) is plotted with the solid line. The imaginary parts of both results are identical as expected, but, 
more interestingly, at the energies of the A* of around 1700 MeV, also the real parts coincide closely. The dashed line 
shows the gauge invariant function G 1 ^ from Eq. (12). This contribution comes from the baryon pole diagrams. As 
there are two baryon propagators, the diagram should be smaller which is indeed the case as Fig. 3 shows. However, 
results can be affected noticeable and one should include this term in general. 

From the comparison in Fig. 3 we see that the cut-off scheme as it has been used in [17] (Kroll-Rudcrman plus 
meson pole term) indeed takes into account the dominant contributions. The baryon pole term, which has not been 
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FIG. 3: Real and imaginary parts of Ai loop functions. Solid line: gauge invariant Gg ; (meson pole term plus Kroll-Ruderman 
term). Dashed line: gauge invariant Ggj (baryon pole term). Dotted line: meson pole plus Kroll-Ruderman term in a cut-off 
scheme (G (cut) + G (cut) ) with A = 881 MeV. The cut-off scheme and Gg ; have identical imaginary parts 



considered in [17], is small. If we would take this term into account in the cut-off scheme, the expression would be 
finite and in the limit A — ► oo, both cut-off scheme and gauge invariant scheme would give identical results. 

Finally, there arc photon couplings to the external baryon of the rescattcring scheme, to vertices of the rcscattcring 
scheme itself, and to intermediate loops of the rcscattcring scheme which all have been ignored in the present study. 
This is because with these couplings, the first loop has no photon attached any longer and is effectively suppressed 
because the integration over the momentum vanishes due to the presence of one s-wave and one p-wave vertex, or 
one s-wave and one d-wave vertex. A detailed discussion can be found in Ref. [11]. There is also the magnetic part 
of the ^B*B* vertex which is proportional to Sa X k [16] with k the photon momentum. As evaluated in Eq. (11) of 
[11], this contribution vanishes for large baryon masses and in practice can be neglected. 



D. Photon coupling to the ttN loop in d-wave 

Diagram (10) of Fig. 2 shows the phototransition amplitude via the irN state in d-wave. The evaluation follows 
the same steps as in Sec. IIB. For irN in (J, I3) = (3/2,1/2), which is the configuration chosen here, |iWr) = 
— y/T/3 \nn + ) + y/2/3 |p7r°). Using standard Feynman rules [17], we obtain 

/ ■+ \(7=3/2,7 3 =l/2) V2 f D + F ~, 

The d-wave meson pole loop function G' Nn has been calculated in Refs. [17] and [29] and we use the results from 
there. In the notation from [17], 

A, ■ f A q 2 1 1 M ( q 2 \ 

Nn J (2tt) 4 (q - kY - m% +te qi- m\ + ie P° - g° - E N (q) + te E N (q) \Ql N (M A .) ) 1 ' 

where QnN{M&*) is the on-shell three-momentum of pion and nucleon for the A* (1700) decay at rest; M is the 
nucleon mass and the other quantities are defined as in Eq. (23). 

The development of a gauge invariant scheme as for the s-wave loops in Sec. IIB, II C is beyond the scope of this 
work. Instead, we follow the lines of Ref. [17]. In the numerical comparison at the end of Sec. II C we have seen 
that the violation of gauge invariance from the cut-off scheme leads almost to the same results as the gauge invariant 
calculation, at least sufficiently beyond threshold; thus, the loop function in Eq. (25) is regularized with a cut-off. 

The cut-off is determined from a comparison between the d-wave meson-baryon loop in dimensional regularization 
with on-shell factorization of the vertices (i.e., as it appears in the rescattering scheme of Sec. II A), and the <i-wave 
loop with cut-off according to 

H^a^-J*^^^ 1 (26) 
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FIG. 4: Effective resonance representation of the radiative decay. 

at the energy y^i of the real part of the resonance position given in Tab. II. In Eq. (26), G„n is the loop in dimensional 
regularization from Eq. (2) with the subtraction constants a^ 7T N) d from Tab. I; Q* N is the on-shell cm. momentum 
from the two e?-wave vertices in on-shell factorization, lu (E) is the pion (nucleon) energy, and M the nucleon mass. 
The cut-off determined in this way is used for the regularization of the meson pole term of Eq. (25). 

E. Effective photon coupling 

In the last sections the amplitudes for the process jp A Hlf °- ) Air have been determined and are written in terms of 
the Tf' 1 ), the unitary solution of the BSE (1) for meson-baryon scattering with the transitions from channel i (A7r, 
H*K, (irN)d) to the Air final state (channel no. 1). In order to determine the partial photon decay width of the 
A*(1700), the T^ 1 ) are expanded around the simple pole in the complex plane with the leading term of the Laurent 
series given by 

\/~S ~ A^A*(1700) 

where gi is the A* (1700) coupling to channel i, the product gig An is provided by the residue and M^*( 1700 ) is the 
complex pole position given in Tab. II. With this replacement for the T^- ll \ the amplitudes from Eqs. (14, 24) can 
be matched to the resonant process shown in Fig. 4, which is given by 

(-it • e) = (-ig A n) r \ f 5A* 7P S 1 * ■ e. (28) 
Vs — Ma* 

This identification allows to write the effective A*(1700)7p coupling (p=proton), gA*-yp, in terms of the one-loop 
photoproduction processes discussed in the last sections: 



5A* 7 p - e SAtt -jj (<^g.i. + <^g.i.J 

C£*K) _ C^*^ NR , j jj s 

5a« 7P - e 9z*k ^ j- (kg.i. + G g i ) (29) 

with the couplings g of the A* (1700) given in Tab. II and the G from Eq. (22). In the same way, the d-wave amplitude 
in Eq. (24) is matched, resulting in 



(jvtt) \/2 D + F ~ 
5a* 7 p = -9{-nN) d e— j — G N7r . (30) 

The effective photon coupling is given by the coherent sum of all processes from Fig. 2, 



5A* 7P - 3a* 7P + 5a* 7P +5a*7 P - \ 6L ) 



III. NUMERICAL RESULTS 

The results for the radiative decay width, given by 



^A'~>yN = T7~ lffA* 7 p| (32) 



k M p 
3tt Ma* 

with <7a*7p from Eq. (31) and k the cm. momentum of the photon, are summarized in Tab. IV. For the 
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TABLE IV: Radiative decay width T of the A* (1700), to be compared with T = 570 ± 254 keV from the PDB [10]. Also, the 
effective couplings g^»^ p of the photon to the A*(1700) via A-7T, H*K, and (ivN)d loops are displayed, in order to show the 
interference pattern of the channels (multiplied with 10 3 ). 





(Air) 




(JVtt) 


r [keV] 


Fit 1 


-82 -i71 


-30-il6 


-9 + i36 


602 


Fit 2 


-85 — z 32 


-28 - i 18 


6 + i34 


403 


Fit 3 


-89-i37 


-25 - i 16 


-1 + i40 


459 


Fit 3' 


-89-i37 


-25-il6 


3 - i41 


730 



experimental value of T = 570 ± 254 keV, we have summed in quadrature the errors from the A* (1700) width and 
the branching ratio into given in the PDB [10]. Table IV also shows the phases of the different contributions. 
The effective couplings from the Att and the ~E*K channels are almost in phase, i.e., they point in a similar direction, 
which makes them dominant over the (ttN)^ channel. 

From the four different results in Tab. IV we prefer the decay width from fit 1. In the other fits some of the free 
parameters have been fixed. As argued in Sec. II A, the remaining free parameters of these extra fits have to absorb 
effects from this reduction of degrees of freedom, and results may become distorted. For fit 3 we have found a very 
similar solution, called fit 3', which has a reversed sign for the coupling of the (7r7V)d-channel to the A* (1700). This is 
a consequence of the weakness of this coupling as explained in Sec. II A. Despite the fact, that under the conditions 
of fit 1, we could not find such an alternative minimum in x 2 , the two different decay widths of 459 keV and 730 keV 
from fit 3 and 3' give an idea of the intrinsic theoretical uncertainties of the present study. Thus, we assign a final 
value of 

T = 602 ± 140 keV 

to the radiative decay width of the A* (1700), with the error given by the difference between fit 3 and 3'. 

The contribution from the T,*K channel is smaller than from the A-7T channel as Tab. IV shows. At first sight this 
seems surprising as g^*K from Tab. II is four times larger than g^. However, the threshold for this channel is at 
y/s = 1880 MeV; at y/s = 1.7 GcV it is closed and the Y,*K loops from Fig. 2 are small. Obviously, the A-7T channel 
with a threshold much lower than the A* (1700) mass dominates the decay; even (Nir) contributes despite to the weak 
coupling to the A*(1700). A similar pattern has been observed in [17] for the radiative decay of the A*(1520) . As a 
consequence, the contribution from the pN channel, which has not been considered in this study, will only moderately 
change the results. 

In the present scheme, the large contributions from the A ++ 7r~ and A°7r + channel add up in the isospin combination 
from Eq. (14), and as a consequence the A-7T channel gives a large contribution to the radiative decay, in good agreement 
with experiment. This is in analogy to the decay A* (1520) — > 7S studied in [17] where the dominant 7rS and 7rE* 
channels add up and result in good agreement with data. 

In Ref. [17], it has also been observed that these channels cancel exactly for the A* (1520) — > 7A decay and the 
discrepancy between experiment and coupled channel model has been attributed to a genuine three quark component 
in the wave function of the A* (1520) , on top of the meson-baryon component from the coupled channel model. For 
the A* (1700) there is no such additional channel for which we can test the model. However, there is a long history of 
calculating radiative decays of excited baryons in quark models, e.g. [30-33]. For the radiative decay of the A* (1700) 
all these works obtain good agreement with experiment. 

Thus, the radiative decay appears well reproduced in both the quark model picture and the present scheme where 
the degrees of freedom are the mesons and baryons. A calculation of electroproduction within the present framework 
could help bring further insight into the question of which theoretical framework is more appropriate; e.g., the quark 
model from [33] has some difficulties in the -D33 channel, but improving the experimental data situation is certainly 
desirable. 



IV. CONCLUSIONS 

In the study of the A* (1700) radiative decay, a model has been formulated in which the photon couples to the 
final loops of the rescattering series that dynamically generates the A* (1700). The irN channel in d-wave has been 
included in the coupled channel scheme. Furthermore, the phototransitions for the dominant s-wave loops have been 
treated in a fully gauge invariant way. 



14 



Previous studies of eleven different pion- and photon-induced reactions have accumulated evidence of a strong 
coupling of the A* (1700) to Arj and E*A', which is a prediction of the unitary coupled channel model. The present 
study provides an extra independent test for the nature of the A* (1700), resulting in additional consistency for the 
picture in which the A* (1700) is dynamically generated. 
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